Quantum Hall Skyrmions with Higher Topological Charge 
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We have investigated quantum Hall skyrmions at filling 
factor v = 1 carrying more than one unit of topological, and 
hence electric, charge. Using a combination of analytic and 
numerical methods we find the counterintuitive result that 
when the Zeeman energy is tuned to values much smaller than 
the interaction energy (gpsB /(e 2 jet) < 9 x 10~ 5 ), the cre- 
ation energy of a charge two skyrmion becomes less than twice 
the creation energy of a charge one skyrmion, i.e. skyrmions 
bind in pairs. The doubly charged skyrmions are stable to 
further accretion of charge and exhibit a 10 % larger spin 
per unit charge than charge one skyrmions which would, in 
principle, signal this pairing. 



Ferromagnetic quantum Hall (QH) systems have quasi- 
particle excitations, "skyrmions", that involve a textur- 
ing of the spins §. There is now strong experimental 
evidence that the lowest energy charged excitations at 
filling factor v = 1 are skyrmions 

A smooth texturing of the spins can be described 
by an effective non-linear sigma model (see ([!]) below) 
where the spin is represented by a unit vector n(r). 
Skyrmions are topological excitations characterized by 
an integer topological charge Z — J drq(r), where q = 
n- (d x nx d y n) / Air is the topological (Pontryagin) density; 
Z is the winding number of the mapping n(r) from the 
compactified plane (S2) to the target space of the sigma 
model (S2 )• A characteristic feature of QH skyrmions is 
that the topological density is proportional to the devia- 
tion of the electron density, p, from its uniform value, p: 
q = v{p — p). Thus, at v = 1, skyrmions with topological 
charge Z carry electric charge —Z\e\. 

Up till now only skyrmions with charge one have been 
considered for the natural reason that one would expect 
Coulombically repelling skyrmions to disfavor charge ag- 
gregation. In this note we study skyrmions at v — 1 with 
higher charge and find that this expectation is violated 
for extremely small values of the Zeeman energy where 
the skyrmions are large objects. More precisely, we find 
that a charge two skyrmion has lower energy than a pair 
of charge one skyrmions for g < g c = 9 x 10~ 5 , where 
g = g^sB I '(e 2 / 'e£) is the dimensionless Zeeman energy. 
This implies that the interaction between two skyrmions, 
which must be repulsive at long distances on account of 
their electric charges, has an attractive core in a region 
where their identities are no longer distinct. We also find 
that skyrmions with charge three and higher are energet- 
ically disfavored at all values of g. 



We begin with variational estimates of the skyrmion 
energies for various topological charges at small g that 
were the motivation for this work. We then study the 
skyrmions numerically, both in the effective sigma model 
by integrating the equations of motion || , and more mi- 
croscopically using a Hartree-Fock scheme jrj . These es- 
tablish the result quoted above and determines g c . Fi- 
nally, we discuss the prospects for observing this effect in 
experiments. 

Effective Sigma Model: The long wavelength physics 
of the spin degrees of freedom in ferromagnetic QH states 
is described by the effective Lagrangian density H 
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Here A is the vector potential of a unit magnetic 
monopole, p s is the spin stiffness (p s = e 2 /(16v / 27re^) 
for v = 1), e is the dielectric constant and I the magnetic 
length. 

The first two terms in C e f / are the leading terms in the 
effective action for any ferromagnet and on account of the 
scale invariance of the gradient term lead, prima facie, 
to charge Z skyrmions of arbitrary size A and energy 
E Z = 4n\Z\p s §. 

The Zeeman and Coulomb terms in £ e ff break the 
scale invariance and their competition sets the size and 
energy of the skyrmions and modifies the detailed form 
of their profiles; these now depend on the dimensionless 
ratio g = (gfisB)/ \e 2 /e£ 2 ) of the Zeeman energy to the 
Coulomb energy. A rotationally symmetric ansatz for a 
skyrmion with topological charge Z is 



VWf) cos(Z9) 



n y = - P{r) sin(Z0) , n z = f(r) 



(2) 



where / obeys the boundary conditions /(0) = 
— 1, /(oo) = 1. The topological density is q — 
[Z /A~Kr)df I dr. The skyrmions in the scale invariant 
sigma model are given by /(r) = \{r/X) 2Z — 4]/[(r / X) 2Z + 
4]. For Z > 1, their topological density has a hollow core. 

Substituting (||) in the action (0) and minimizing leads 
to a non-linear, non-local integro-differential equation for 
/. Below, we discuss results obtained by numerically in- 
tegrating this equation using a relaxational technique H . 
But first, it is instructive to consider approximate solu- 
tions at small g; the details are discussed elsewhere [gj. 
For Z > 1 we take the solution to be of the form of the 
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scale invariant sigma model solution with an optimized 
scale parameter A. This yields, 

A = 0.780ftT 1/3 and E = ^ [2y^ + 2.86g 1/3 ] [Z = 2) 
A = 1.29£g~ 1/3 and E = ^ [3y^ + 4.82p 1/3 ] (Z = 3) 
A = 1.75^- 1/8 and £ - ^ [4^ + 7.2I5 1 / 3 ] (Z = 4) 

for the lowest three values of Z. Note that the energy per 
unit charge is the same for the g independent piece as is 
appropriate for solutions of the scale invariant problem 
]S| while for the g dependent piece it increases monoton- 
ically with Z. Hence, there is no binding between this 
set of skyrmions. 

However, this procedure runs into trouble with the 
Z = 1 skyrmion where the scale invariant solution yields 
a Zeeman energy that diverges logarithmically with sys- 
tem size for any A. This can be fixed by matching the 
scale invariant solution onto the exact asymptotic solu- 
tion in the outer region which decays exponentially with 
a length oc lg~ x / 2 . As a consequence we get logarthmi- 
cally modified expressions pL 
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= 1, c mcr |0) = and & m(J creates 



A = 0.558^(g|lng|)- 1/3 and 

^=^[^ + 0.622(.g|ln3|) 1/3 ] (Z = 1) , (4) 

for the core scale parameter and energy. It is evident 
that the presence of the logarithm by itself implies that 
charge one skyrmions will pair bind for small enough g 
and that the binding energy will vanish as g vanishes. 
Using the above expressions one finds a critical value of 
g c w 5.3 x 10~ 6 for this to happen. That such a critical 
value has to exist can be inferred by a separate micro- 
scopic computation of the pair-binding energy at large g, 
where the skyrmion size is 0(£) and they reduce to polar- 
ized quasiparticlcs, which shows the abscence of binding. 
However, the quantitative issue is rather delicate. Even 
at small g the result (^) is expected to have logarithmi- 
cally subdominant corrections and is known numerically 
to be not terribly accurate in estimating the difference 
E(g) — E(0), even for g as small as 10~ 6 M. Also, the 
asymptotic validity of our small g expressions is not on 
completely rigorous footing at present. So it follows that 
while the pair binding should be present only at small 
<7, a reliable estimate of the critical g c requires the more 
accurate computations of the skyrmion energies reported 
below. 

Hartree-Fock: We complement the effective action cal- 
culations by the Hartree-Fock method introduced by Fer- 
tig et. al. 0]. For skyrmions with topological charge Z 
the Hartree-Fock wave function is fixed by requiring it to 
be in the lowest Landau level and by the symmetries of 
the classical solution (@) to be 



m=-Z 

where \u m \ 2 + \i 
electrons with spin a in the angular momentum m-state 
)in the lowest Landau level. We choose u m , v m real and 
require u m —* 1 as m — > 00 . The coefficients u m , v m 
are determined by numerically iterating the Hartree-Fock 
equations until a self-consistent solution is found. Obvi- 
ously, Y$>) has electric charge —Z\e\ relative to the polar- 
ized groundstate \^>q) = rim=o c mTl^)- Skyrmions with 
charge Z < are often referred to as antiskyrmions. In 
the actual numerical calculation we have a finite num- 
ber of particles, m < N, and impose un — 1 (or 
un = uat_i = 1 see below). 

Results: Our results are shown in Fig. 1-3. In Fig. 1 we 
give the energies to create one Z = 2 skyrmion and to 
create two Z = 1 skyrmions, both at fixed particle num- 
ber and magnetic field, as functions of g. These energies 
are equal for the skyrmion (Z > 0) and the antiskyrmion 
(— Z < 0). This means that we have shifted the Hartree- 
Fock result for the state (||) by a constant. For large g, 
the Z = 2 skyrmion has higher energy than two Z = 1 
skyrmions. In this region the skyrmions are small, much 
smaller than the system size, and the Hartree-Fock cal- 
culation is reliable. 
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FIG. 1. Energies for one Z = 2 skyrmion and for two Z = 1 
skyrmions. The lines are Hartree-Fock results for N = 1000 
and the points are sigma model results. The spin of one Z = 2 
and two Z = 1 skyrmions is indicated along the lines. The 
inset shows that the energies cross at g c — 9 x 10 -5 
model results only). E(0) is the g = energy. 



(sigma 



The sigma model is to be trusted only for small g when 
the skyrmions are large. It nevertheless gives energies 
that are close to, but slightly larger than, the Hartree- 
Fock energies for the range of g shown. As g decreases, 
the difference in energy between one Z — 2 and two 
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Z = 1 skyrmions decreases. We see that, except pos- 
sibly for very small g, the Z = 2 skyrmion has larger en- 
ergy. As noted earlier, in the limit g — > the skyrmions 
become infinitely large and only the sigma model term 
in (|l|) contributes to the energy. Thus both energies 
must eventually approach the pure sigma model result 



E(0) = 87175 s 



/e£ as g — > 0. To investigate the 



region of small g, we proceed in two different ways. 

The numerical sigma model method can be used for 
small g. The inset in Fig. 1 shows the difference in en- 
ergy as a function of g obtained using the sigma model. 
This predicts a crossing at g c = (9 ± 1) X 10~ 5 . At this 
point the energy is Eg c ~ 0.751e 2 /e£ and there is a 10% 
change in the spin per charge (S/Z) of the lowest energy 
excitation: (S/Z) z =i = 221 and (S/Z) z=2 = 243. The 
size of the excitation also changes: The charge radii (per 
unit charge) are rx=i = 21£ and r z =2 — 31£ We can also 
see that when g decreases below g c , the difference in en- 
ergy first increases but then eventually starts to decrease 
as it must, since it should vanish at g = 0. 

Using Hartree-Fock for smaller g than shown in Fig. 
1, we find apparent crossing points, g c (N). However, 
the g c (N) are so small that the skyrmions are large and 
there are large finite size effects: g c depends on N and on 
whether one imposes as boundary condition that one or 
several of the spins are up at the edge of the system. It is 
difficult to increase N very much beyond N = 1000 since 
the convergence is very slow for the small values of g we 
are interested in. We conclude that within Hartree-Fock 
we are not able to directly probe a region of g where an N- 
independent crossing might take place. However, we can 
reach smaller g by finite size scaling. In Fig. 2 we plot 
the crossing point, g c (N), as a function of 1/N for three 
different sets of boundary conditions (as indicated in the 
figure) JTo[ ] . To each set we fit a quadratic polynomial in 
1/N and read off g c = g c (oo) = (8.9 ± 1.4) x 10~ 5 . The 
error is the standard deviation in the three extrapolated 
values. The errors in the fit for each curve are much 
smaller. Also, adding or subtracting a point at small 
N does not affect the results. That the three different 
sets, which give very different g c (N), predict virtually 
the same nonzero g c indicates strongly that there is a 
crossing. This value of g c agrees with the one obtained 
in the effective theory. Thus we conclude that the doubly 
charged skyrmions have lower energy for g < g c = 9 x 
10- 5 . 

It is straightforward to calculate not only the energy 
but also the total spin as well as spin density and charge 
density. The total spin, which is indicated along the line 
in Fig. 1, is much larger for the doubly charged skyrmion 
than for the charge one skyrmion. In Fig. 3 we give 
examples of densities for the doubly charged skyrmion 
and compare to one charge one skyrmion. Note that the 
charge density has its maximum a finite distance away 
from the origin for the doubly charged skyrmion (as the 
scale invariant skyrmion has). 
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FIG. 2. Finite size scaling of g c (N) using Hartree-Fock re- 
sults for antiskyrmions with various boundary conditions (one 
or two spins up at the edge, as indicated in the figure). g c {N) 
is where the energy of one Z — — 2 skyrmion is equal to the 
energy of two Z = — 1 skyrmions in a Hartree-Fock calcu- 
lation with N particles. Note that g c (N) depends strongly 
on N and on the boundary condition. Fitting to quadratic 
polynomials in 1/N gives g c = g c (oo) = 8.9 x 10~ 5 . 
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FIG. 3. Charge and spin densities at g = 0.02, 0.01, 0.005 
for one Z = 2 skyrmion (using Hartree-Fock). Densities for 
one Z = 1 skyrmion at g — 0.02 is included for comparison. 

The numerical sigma model gives lim^ol-^ ~ 
y ^T/8e 2 /e£]/g 1 ^ = 2.86e 2 /e£ for Z = 2 in agreement 
with (E3) . Thus there is no logarithm in the Z — 2 energy 
and hence no reason to expect that Z > 2 skyrmions 
would ever have the lowest energy per charge. We have 
also checked that Z = 3 skyrmions have higher energy 
(per charge) than Z — 1, 2 skyrmions for g > 10~ 7 . 
Discussion and Experimental Implications: Our 
analysis has established that two skyrmions are unsta- 
ble to pair-binding at very small Zeeman energies. At 
distances much larger than their size, the inter-skyrmion 
potential must be repulsive as it will be dominated by 
their Coulomb interaction. It follows that the attraction 
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that we have found sets in only when the skyrmions over- 
lap and their identities begin to merge. In computing the 
skyrmion energies, we have ignored the effects of Landau 
level mixing and of the softening of the Coulomb interac- 
tion at small distances by the finite extent of the subband 
wavefunctions transverse to the plane of the electron gas. 
It is likely that neither of these makes a huge difference to 
the energetics. In the sigma model description, the domi- 
nant effect of Landau level mixing is a renormalization of 
the spin stiffness which does not affect the inter-skyrmion 
energetics. The softening of the Coulomb interactions is 
significant at distances which are 0(£) and for much big- 
ger skyrmions (such as in the region of interest where 
they exceed 20£) most of the self-interaction is unsoft- 
ened. (We also note that the long wavelength arguments 
regarding the peculiarities of the Z = 1 solution that un- 
derlie the pairing depend only on the form of the sigma 
model action; they will apply mutatis mutandis to other 
ferromagnetic fillings as well.) 

This is encouraging in that it suggests that it might be 
possible to see the binding of the skyrmions as a sharp 
decrease in the spin polarization with decreasing g in the 
vicinity of v = 1, where the ground state will exhibit a di- 
lute density of skyrmions. Unfortunately, this is unlikely 
on two grounds. First, it would require an extremely 
fine tuning of g which seems hard to accomplish. Sec- 
ond, there is the more fundamental problem that when 
the skyrmions get to be very big, they become sensitive 
to the details of the disorder potential which will tend 
to limit their size. Effectively, we expect the disorder to 
provide a lower limit to g and while we have not com- 
puted this for realistic disorder realizations, we are not 
greatly encouraged by the data in Q which found a max- 
imum skyrmion size, even for g reduced under pressure, 
considerably smaller than the critical one calculated in 
this paper. 

Note Added: While finishing this work we became 
aware of a recent paper by Nazarov and Khaetskii fTl| ] 
that studies higher topological charge skyrmions in the 
effective sigma model in the spirit of our variational cal- 
culations. These authors also argue that the interaction 
between two charge one skyrmions is attractive at inter- 
mediate and short distances and that the binding dis- 
tance is zero, at which they form a charge two skyrmion. 
However, their quoted results are independent of g, in 
contradiction with ours which hold only at small g, and 
they appear not to have noticed the existence or the ex- 
treme smallness of g c . 
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